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SECTION  I 


INTRODUCTION 

The  Longley-Rice (  1 *(20)  and  j0hnson-Gierhart(t<)“(t) 
prediction  programs  of  the  Institute  for  Telecommunication  Sciences 
nre  statisticnl/semi-empir ica]  models  of  tropospheric  radio 
propagation  for  low  snd  high  altitude  scenarios  respectively  in  the 
frequency  range  0.0?  -  i»0  GHz.  These  programs  are  restricted  to 
frequencies  above  ?P  MHz  because  "sky"  and  "ground"  wave  propagation 
paths,  which  can  be  dominant  propagation  paths  at  frequencies  less 
than  ?P  MHz,  are  not  included  in  these  programs.  These  programs  are 
restricted  to  frequencies  less  than  R0  GHz  because  the  empirical  data 
base  does  not  include  absorption  and  refractivity  of  the  atmsophere  or 
ground  at  wavelengths  shorter  than  1  cm. 

Version  1.?.1  of  the  Longley-Rice  program  and  the  Air-to-Air 
(ATOA)  version  of  the  Johnson-Gierhart  program  have  been  acquired  by 
mttre  Corporation  for  prediction  of  propagation  path  loss  in  radio 
scenarios.  Version  1.2.1  of  the  Longley-Rice  program  is  written  in 
lOffi  ANSI  Fortran  and  therefore  is  compatible  with  most  large 
computers.  The  Johnson-Gierhart  ATOA  program,  which  was  originally 
written  for  a  CDC  CYBER-170/750  computer  has  been  converted  and  made 
compatible  with  the  IBM-STO/^O’I  computer  at  MITRE  Corporation. 

The  theory,  computer  programs,  and  user's  guides  for  the 
Longley-Rice  and  Johnson-Gierhart  prediction  models  are  given  in 
References  Ml  -  f 6 1 ,  F201.  This  paper  discusses  the  input  and  output 
parameters  of  these  prediction  models,  with  particular  emphasis  on 
input  parameter  specification. 
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The  outputs  or  these  programs  are  statistical  values  of  nbasic 
transmission  loss".  In  Section  2,  basic  transmission  loss  is  defined 
and  the  particular  quantities  of  baslo  transmission  loss  which  are 
evaluated  by  these  programs  are  discussed  and  oompared  with 
theoretical  values.  The  motivation  for  using  these  programs  are  also 
given. 

Input  parameter  specifications  for  the  Longley-Rice  and 
Johnson-Glerhart  prediction  models  are  described  in  Sections  3  and  4 
respectively. 
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SECTION  2 


SASIC  TRANSMISSION  LOSS 


With  reference  to  Figure  1,  the  predetection  system  margin 
M(d,r.)  (in  dB)  of  a  radio,  for  a  great  circle  path  distance  d  and 
message  quality  rj,  is  given  by 

S(d) 

f - ^ - \ 

*  Pt  -  LH.T  '  LH.T  -  LC.T  *  °T  -  V11*  *■>(,-"-  Br<rl>  <*•'> 

where 


S(d)  s  available  signal  power  at  the  output  terminals  of  the 
equivalent  lossless  receiving  antenna  (dBm) 

s  transmitter  carrier  available  outut  power  (dBm) 

L„  _  =  insertion  loss  of  transmitter  transmission  line 
’  (including  reflection  losses,  if  any)  (dB) 


ohmic  loss  of  matching  network  for  the  transmitting 
antenna  (dB) 


lc  T  =  ohmic  loss  of  the  transmitting  antenna  (dB) 

Dt  s  directivity  of  the  transmitting  antenna  (dBi) 

Lfe(d)  s  basic  transmission  loss  of  propagation  path  (dB) 

Dr  s  directivity  of  the  receiving  antenna  (dBi) 

N  :  system  available  noise  power  at  the  output  terminals  of 
the  equivalent  lossless  receiving  antenna  (dBm) 

R  (r, )  :  required  predetection  slgnal-to-nolse  ratio  (dB) 

r  x 
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FIGURE  1:  RADIO  SYSTEM  MARGIN  PARAMETERS 


The  convention  is  followed  that  capitalized  margin  parameters  are  in 
units  of  dB  whereas  uncaptiallzed  margin  parameters  refer  to  their 
numerical  values. 

The  basic  transmission  loss  L^(d)  is  the  margin  parameter  random 
variable  in  Eq.  (?.1)  that  is  a  function  of  the  propagation  range  of 
great  circle  distance  d  between  the  transmitting  and  receiving 
antennas.  More  specifically,  "the  basic  transmission  loss  (sometimes 
called  path  loss)  of  a  radio  circuit  is  the  transmission  loss  expected 
between  ideal,  loss-free,  isotropic,  transmitting  and  receiving 
antennas  at  the  same  locations  as  the  actual  transmitting  and 
receiving  antennas"'  . 

The  basic  transmission  loss  L^(d)  (in  dB)  may  be  expressed  as: 

4)(d)  *  So*'0  *  U(d)  *  V(d)  (2.2) 

where 

L^0(d)  *  local  path  propagation  loss  for  a  path  with  no  buildings 
or  significant  vegetation  in  the  immediate  vicinity  of 
the  antennas  (hourly  median  value,  in  dB) 

U(d)  :  urban  area  propagation  loss  term  resulting  from 

buildings  in  the  Immediate  vicinity  of  the  antennas 
(hourly  median  value,  in  dB) 

V(d)  s  vegetative  propagation  loss  term  resulting  from 

significant  vegetation  in  the  immediate  vicinity  of  the 
antennas  (hourly  median  value,  in  dB) 

The  propagation  path  loss  given  by  Eq.  (2.2)  does  not  Include 
rapid  fading,  about  the  hourly  median  value,  of  the  received  signal 
which  results  from  multipath  interference  (vertical  loblng)  for  a 
small  change  in  range  (of  the  order  of  a  oarrier  wavelength)  in  a 
mobile  system  operating  over  Irregular  terrain.  Instead  of 
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superimposing  the  distributions  of  hourly  median  and  rapid  fading 

propagation  losses,  the  effeots  of  multipath  rapid  fading  oan  be 

included  in  the  system  margin  model  by  increasing  the  value  of  the 

required  predeteotion  signal-to-noise  ratio  R  (r.)  necessary  to  obtain 

*  1  (8) 

a  specified  channel  quality  r^  in  the  presenoe  of  rapid  fading 

The  local  path  propagation  loss  Lfeo(d)  in  Eq.  (2.2)  may  be 
expressed  as: 

Lbo(d)  a  Lbf(d)  ♦  A(d)  (2.3) 

where 

L.  _(d)  s  free-spaoe  propagation  loss  (dB) 


=  32.W  ♦  20  log1()  fMHz  ♦  20  log10  d^ 


d^  =  great  circle  distanoe  between  transmitter  and  receiver 
antennas  (km) 

X  a  RF  carrier  wavelength  (in  units  of  d) 

fMHz  s  carrier  frequency  (MHz) 

c  *  free-space  velocity  of  propagation  =  0.29979  km/MS 

A(d)  a  excess  propagation  loss  over  that  of  free  space  for  a 

path  with  no  buildings  or  significant  vegetation  in  the 
immediate  vicinity  of  the  antennas  (hourly  median  value 
in  dB).  This  term  is  usually  modelled  by  semi-empirical 
methods . 

The  Lor.gley-Rloe  and  Johnson-Gierhart  prediction  programs  are 
concerned  with  estimating  the  hourly  median  values  of  excess 
propagation  loss  A(d)  defined  by  Eq.  (2.3).  The  urban  loss  U(d)  and 
vegetative  loss  V(d)  terms  in  Eq.  (2.2)  are  not  predicted  by  these 
programs.  The  rapid  fading,  about  the  hourly  median  values,  is  also 
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not  predicted  by  the  Longley-Rice  program  but  Is  an  available  vertical 
loblng  option  In  the  Johnson-Gierhart  program.  Before  discussing  the 
excess  loss  A(d),  a  brief  review  is  given  here  of  the  urban  loss  U(d) 
and  vegetative  loss  V(d)  terms  in  Eq.  (2.2). 

The  urban  area  propagation  loss  term  U(d)  in  Eq.  (2.2)  nas  been 

(9) 

estimated  to  have  a  median  value  U(d,50f)  given  by  . 

16.5  ♦  15  lo810  ”  0,12  dkm’  urDan  area 

U(d,50%)  =  (2.9) 

0,  both  antennas  are  in  open  areas 

which  agrees  within  1  dB  with  empirical  data  at  100  -  300  MHz  and 
distances  10  -  70  km.  The  median  value  U(d,50f)  given  by  Eq.  (2.4)  is 
the  difference  in  the  median  values  of  excess  propagation  loss 
reported  by  Okumura(10^  for  urban  areas  and  by  Longley-Rlcet2)  for 
open  areas.  Eq.  (2.4)  is  based  on  data  for  which  the  receiver  antenna 
was  near  ground  level  (at  a  height  of  3m)  and  the  transmitter  antenna 
was  at  various  elevated  heights  of  30  -  600  m.  At  a  frequency  of  88 
MHz  and  a  distance  of  35  km,  the  median  value  of  the  additional 
transmission  loss  from  urban  area  clutter  is  found  from  Eq.  (2.4)  to 
be  11.5  dB, 

The  vegetative  propagation  loss  term  V(d)  in  Eq.  (2.2)  is 
appreciably  less  at  VHF  frequencies  than  at  higher  frequencies  because 
vegetation  is  appreciably  more  transparent  at  longer  wavelengths  and 

because  obstacles,  suoh  as  vegetation,  diffract  more  energy  into 

(0) 

shadow  zones  at  longer  wavelengths  .  Only  vegetation  in  the 
immediate  vicinity  of  the  antennas  should  be  considered  in  estimating 
V(d)  in  Eq.  (2.2)  because  knife-edge  diffraction  by  vegetation  distant 
from  the  antennas  is  usually  included  in  the  seml-empirioal  methods 
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used  for  estimating  the  excess  propagation  loss  A(d).  The  loss  term 
V(d)  Is  the  lessor  of  the  absorptive  path  loss  through  the  vegetation 
and  the  diffractive  path  loss  over  the  vegetation. 

When  one  of  the  antennas  is  placed  near  a  grove  of  trees  or  in  a 

Jungle,  vertically  polarized  VHF  radio  waves  are  attenuated 

appreciably  more  than  horizontally  polarized  waves.  For  example,  at 

30  MHz  and  100  MHz,  the  average  loss  from  nearby  trees  was  reported  to 

be  2  -  3  dB  and  5  -  10  dB  respectively  with  vertical  polarization  and 

approximately  0  dB  and  2  -  3  dB  respectively  for  horizontally 

polarized  signals'1^  *  ^ 1 1  \  In  dense  jungles,  vertically  polarized 

waves  can  be  attenuated  about  15  dB  more  than  horizontally  polarized 
(o)  (1?) 

fields  '  .  At  higher  frequencies,  the  effects  of  polarization  on 

vegetative  loss  are  not  as  pronounced. 

In  experimental  studies  by  La  Grone^1^  of  propagation  of 
horizontally  polarized  waves  behind  a  grove  of  3  m  tall  live-oak  and 
blackberry  trees  on  flat  ground  in  Texas  at  frequencies  0.5  -  3  GHz 
end  at  distances  greater  than  five  times  the  tree  height,  measurements 
of  path  loss  were  in  good  agreement  with  theoretical  oredictlons  of 
diffraction  over  an  ideal  knife  edge  assuming  distances  and  heights 
the  same  as  those  in  the  measurements.  For  such  a  case,  the  loss  term 
V(d)  may  be  interpreted  as  the  difference  in  losses  between  knife-edge 
(tree)  diffraction  and  smooth  spherical  earth  diffrsotion  with  losses 
expressed  in  dB.  Approximate  numerical  values,  deduced  from  data  for 
the  above  case,  are  V(d)  =  -H,  -2,  and  *2  dB  for  reoelver  heights 
above  local  terrain  of  2,  10,  and  18  m  respectively,  e  frequency  of  82 
MHz,  transmitter  to  reoelver  distanoe  d  *  67  km,  transmitter  height  of 
U2H  m,  and  receiver  to  grove  distance  of  111  m. 
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Estimates  of  the  excess  propagation  loss  A(d)  defined  in  Eq.  (2.3) 
can  be  as  formidable  to  calculate  as  the  estimates  for  urban  area  and 
vegetative  propagation  losses  because  of  the  serai-empirical  nature  of 
the  required  models. 

The  excess  propagation  loss  A(d)  Is  generally  a  stochastic 
quantity  because  the  scenarios  of  interest  are  generally  not  for 
deterministic  propagation  paths  but  are  for  specified  classes  of 
propagation  paths.  For  example,  the  propagation  paths  may  be 
specified  as  being  over  irregular  terrain  characterized  by  rolling 
plains,  average  ground  permittivity,  and  random  siting. 

Radio  waves  generally  may  be  propagated  (a)  through  or  along  the 
surface  of  the  earth  (ground  wave),  (b)  through  the  lower  atmosphere 
of  the  earth  beneath  the  ionosphere  (tropospheric  propagation)  or  (c) 
by  reflection  or  scatter  In  the  upper  atmosphere  (sky  wave)  from 
natural  reflectors  (Ionosphere,  aurora)  or  artificial  reflectors 
(satellites).  At  frequencies  greater  than  20  MHz,  ground  wave 
propagation  losses  (except  for  very  short  paths  within  the  radio 
horizon  and  along  the  earth's  surface)  and  sky  wave  propagation  losses 
(except  for  very  long  propagation  paths  beyond  the  radio  horizon)  are 
usually  very  muoh  larger  than  tropospheric  propagation  losses.  The 
Longley-Rice  and  Johnson-Gierhart  programs  consider  only  tropospheric 
propagation  paths. 

For  tropospheric  propagation  over  irregular  terrain,  the  possible 
modes  of  propagation  may  be  oategorlzed  as: 

(1)  multipath  Interference 

(2)  multipath  -  diffraction  transition 

(3)  diffraction  (smooth  spherical  earth  and  knife-edge) 

(U)  diffraction  -  tropospheric  soatter  transition 

(5)  tropospherle  soatter  . 


The  regions  for  these  propagation  inodes  are  shown  in  Figure  2. 

Node  (1)  is  the  dominant  mode  of  propagation  for  line-of-sight 
paths  which  clear  the  radio  horizon  by  greater  than  approximately  1/8 
of  a  Fresnel  number  where  the  Fresnel  number  is  the  number  of 

half-wavelengths  of  the  path  difference  between  the  direct  ray  and  the 

(IH) 

indirect  ray  which  is  specularly  reflected  from  the  ground  .  Mode  (2) 
occurs  for  line-of-sight  propagation  paths  which  are  within  1/8 
Fresnel  number  of  the  radio  horizon.  Mode  (3)  occurs  for  propagation 
paths  which  are  beyond  the  radio  horizon  by  more  than  1/8  Fresnel 
number  but  less  than  that  for  which  tropospheric  scatter  starts  to 
become  significant.  Mode  (H)  is  a  transition  mode  between  diffraction 
and  troposeatter  modes.  Mode  (5)  occurs  for  propagation  paths  which 
are  sufficiently  beyond  the  radio  horizon  that  tropospheric  scatter 
losses  are  less  than  diffractive  losses.  Except  for  mode  (1)  lobing, 
the  excess  propagation  loss  A(d)  generally  Increases  with  decreasing 
height  h^  as  the  dominant  mode  of  propagation  progresses  from  (1)  to 
(5). 

As  an  example  of  tropospheric  modes  (1)  and  (2),  consider  the 
various  scenarios  shown  in  Table  1.  For  smooth  earth,  all  of  the 
Table  1  scenarios  correspond  to  radio  links  within  or  on  the  radio 
horizon.  However,  for  random  siting  of  the  ground-based  radio  on 
Irregular  terrain,  the  radio  line-of-sight  to  low  altitude  aircraft 
will  often  be  obstructed.  In  the  case  of  nonobstructed  radio 
line-of-sight  over  smooth  terrain,  the  principal  mode  of  propagation 
is  smooth  spherical  earth  diffraction  coupled  with  multipath 
Interference  between  the  direct  and  indirect  signals  reflected  by  the 
terrain  to  the  receiver.  In  the  case  of  an  obstructed  radio 
line-of-sight  over  Irregular  terrain,  the  principal  mode  of 
propagation  is  smooth  spherical  earth  dlffraotlon  ooupled  with 
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FIGURE  2:  REGIONS  FOR  VARIOUS  TROPOSPHERIC  PROPAGATION  MODES  OVER  IRREGULAR  TERRAIN 


HASH  TRAM-MISSION  LOSS.  U  'd®) 


TABLE  1:  BASIC  TRANSMISSION  LOSS  FOR  LOW  AND  HIGH  ALTITUDE  SCENARIOS,  f  =  50  MHz 


knife-edge  diffraction  by  the  obstructing  terrain.  Multipath 
interference  is  the  dominant  mode  of  propagation  in  scenario  numbers 
3b,  7b1,  and  8  which  dear  the  radio  horizon  for  a  smooth  earth  by  at 
least  1/8  Fresnel  number.  All  of  the  other  scenarios  in  Table  1, 
which  clear  the  radio  horizon  by  less  than  0.01  Fresnel  number, 
correspond  to  the  multipath  -  diffraction  transition  mode. 

In  the  multipath  -  diffraction  transition  mode,  the  propagation 
path  loss  is  significantly  larger  than  the  free  space  loss.  In  the 
case  of  multipath  interference,  the  interference  is  almost  totally 
destructive  because:  1)  the  Fresnel  amplitude  reflection  coefficient 
is  approximately  -1  at  low  grazing  angles  of  incidence;  2)  the  surface 
roughness  reflection  coefficient  is  approximately  unity  at 
sufficiently  low  grazing  angles  of  incidence;  3)  the  path  length 
difference  between  the  direct  and  indirect  signals  Is  much  less  than  a 
wavelength . 

In  the  case  of  diffraction,  the  path  loss  increases  exponentially 
with  increasing  distance  of  the  transmitter  or  receiver  into  the 
shadow  region  of  the  obstructing  terrain.  For  a  given  distance  in  the 
shadow  region,  the  sharpness  of  the  obstructing  terrain  appreciably 
alters  the  path  loss  in  a  diffraction  mode.  Therefore,  for  a 
diffraction  mode  of  propagation,  slope  and  height  distribution  of  the 
obstructing  terrain  affect  the  path  loss. 

The  Longley-Rloe  and  Johnson-Gierhart  statistical,  seml-emplrlcal 
programs  are  particularly  useful  in  modeling  propagation  loss  over 
irregular  terrain  in  the  transition  modes  (2)  and  (4)  of  Figure  2. 

The  propagation  loss  for  mode  (2)  is  found  from  empirical  data  snd 
from  extrapolations  between  theoretical  models  for  multipath 
interference  and  smooth  spherloal  earth  diffraotlon.  The  propagation 
loss  for  mode  (<t)  is  found  from  empirical  dats  and  from  extrapolations 
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between  theoretical  models  for  smooth  spherical  earth  diffraction  and 
tropospheric  scatter . 

The  Longley-Riee  prediction  program  is  applicable  to  scenarios 
where  both  the  transmitter  and  receiver  antennas  are  at  heights  above 
local  ground  between  0.5  m  and  3  km  (which  we  shall  designate  as  "low 
altitude"  scenarios').  The  Johnson-Gierhart  program  is  applicable  to 
"high  altitude"  scenarios  in  which  (1)  the  lower  antenna  is  at  a 
height  above  local  ground  between  0.5  u  and  approximately  3  km;  (2) 
the  higher  antenna  is  less  than  100  km  but  at  a  sufficient  height 
ahove  local  ground  that  the  elevation  angle  at  the  lower  antenna  of 
the  terrain  limited  radio  horizon  is  less  than  the  elevation  angle  of 
the  higher  antenna;  and  (3)  the  terrain-limited  radio  horizon  for  the 
higher  antenna  is  taken  either  as  a  common  horizon  with  the  lower 
pntenna  or  as  a  smooth  earth  horizon  with  the  same  elevation  as  the 
lower  antenna  effective  reflecting  plane.  These  altitude  restrictions 
and  the  use  of  these  programs  are  based  on  the  following 
considerations : 

a.  Whereas  two-ray  multipath  interference  models  are  adequate 

for  path  clearances  greater  than  1/8  Fresnel  number  and  whereas  smooth 

earth  spherical  diffraction  models  are  adequate  for  transhorizon  paths 

well  beyond  the  radio  horizon,  an  extrapolation  between  these  models, 

even  for  a  smooth  earth,  is  presently  required  for  modeling  of 

(18) 

propagation  paths  near  the  radio  horizon  ,  Reference  [lit]  gives  a 
deterministic  computer  program  for  such  an  extrapolation.  However, 
for  path  loss  averaged  over  random  paths  above  irregular  terrain  near 
the  radio  horizon,  a  semi -empirical ,  stochastic  extrapolation  is 
required.  The  empirical  weighting  accounts  for  knife-edge  diffraction 
effects  over  a  rough  earth.  The  Longley-Rice  semi-empirical 
prediction  program  does  such  an  extrapolation  and  allows  for  double 
horizon  diffraction  for  both  random  and  specific  terrain  profiles. 
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b.  Probabilistic  predictions  of  path  loss  are  possible  because 
the  data  base  inoludes  many  samples  for  various  locations,  time  of 
year,  and  experimental  situations.  Much  of  the  data  base  for  the 
Longley-Rlce  program  Is  for  double  horizon  diffraction  paths  and  was 
obtained  In  the  frequency  range  20  -  100  MHz*1^.  Much  of  the  data 
base  for  the  Johnson-Glerhart  program  is  from  200  single  horizon 
diffraction  paths  contained  in  the  data  of  Longley,  et  al^\ 

c.  The  Longley-Rlce  program  assumes  a  uniform  atmosphere  (linear 
refractive  gradient)  and  Is  therefore  not  applicable  to  propagation 
paths  in  standard  exponential  atmosphere  at  elevations  above  3000  m. 
The  Longley-Rlce  program  is  also  restricted  at  each  antenna  to  path 
elevation  angles  less  than  12°.  For  path  elevation  angles  greater 
than  12°,  time  variability  of  path  loss  caused  by  atmospheric 
refraction  is  appreciably  less  than  that  of  the  empirical  data  base 
(which  Is  limited  to  refractive  effects  at  elevation  angles  less  than 
12°). 


d.  The  Johnson-Glerhart  prediction  program  is  restricted  to 
single-horizon  diffraction  which  allows  for  ray  tracing  In  standard 
atmospheres  from  the  horizon  back  to  the  antenna  site.  The 
Johnson-Glerhart  program  Is  therefore  applicable  to  paths  at  high 
elevations  and  steep  elevation  angles  but  Is  not  applicable  at  low 
elevations  where  double  horizon  diffraction  may  be  significant. 

In  the  Longley-Rlce  and  Johnson-Glerhart  programs,  ionospheric 
propagation  and  auroral  soatter  effects  are  assumed  to  be  negligible. 
This  assumption  is  valid  at  frequencies  above  100  MHz  and  is  expected 
to  be  valid  at  frequencies  20  -  100  MHz  for  sufficiently  short  path 
lengths.  In  the  Longley-Rlce  program  a  print -out  warning  is  issued  at 
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frequencies  less  tnan  *»0  MHz  to  remind  tne  user  tnat  sxy  wave  enecua 
may  be  important  at  sufficiently  long  patn  lengths,  me 
jonnson-Uiernart  program  can  be  extended  to  frequencies  as  low  as 
MHz  provided  tnat  sky  wave  effects  are  negligible. 

Both  of  tnese  prediction  programs  restrict  antenna  sites  to 
locations  lor  which  the  ratio  of  the  distance  to  the  terrain-limited 
radio  horizon  to  that  lor  a  smooth  spnencai  eartn  is  greater  tnan  0.1 
and  less  than  i.Q.  This  restriction  applies  to  both  antenna  sites  in 
the  Longiey-Rice  program  ano  only  to  tne  lower  antenna  in  tne  AlUA 
program,  tor  example,  these  programs  would  not  be  applicable  to  a 
scenario  in  whicn  tne  ground  site  is  at  the  bottom  of  a  steeply  rising 
hill  because  the  welgnted  extrapolation  of  models  in  these  programs  do 
not  allow  for  severe  knife-edge  diffraction.  Tne  hourly  median  Dasic 
transmission  loss  computed  oy  these  programs  does  not  include  losses 
resulting  from  ionage  or  buildings. 

In  Figure  j,  the  cumulative  distribution  functions  of  path  loss 
are  plotted  on  normal  probability  paper  lor  both  nigh  and  low  altitude 
scenarios,  me  path  loss  tor  the  low  altitude  scenario  of  Table  l 
(scenario  no.  9a)  was  predicted  utilizing  the  Longiey-Rice  model 
whereas  the  path  loss  lor  the  mgn  altitude  scenario  (scenario  no. 

7c2)  was  predicted  utilizing  the  Johnson-Uierhart  model.  Whereas 
Longiey-Rice  program  predictions  of  patn-  loss  for  low  altitude 
scenarios  are  approximately  normally  distributed  when  path  loss  is 
expressed  in  dB,  johnson-uierhart  program  predictions  of  path  loss  tor 
high  elevation  scenarios  are  two-piecewise  normally  distributed  with  a 
breakpoint  at  the  median  value  as  explained  oelow.  Location,  time, 
and  situation  (model)  uncertainties  contribute  to  path  loss 
variability,  in  the  Longiey-mce  program,  location  variability  is 
usually  the  dominant  path  loss  variability  because  diffraction  by 
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JOHNSON -GIERHART  (h.  *  2.5  »,  h,  =  4600  tn,  d  =  225  km) 


terrain  at  both  ends  of  the  link  Is  considered.  In  the  Johnson- 
Gierhart  program,  time  variability  is  usually  the  dominant  path  loss 
variability  because  diffraction  by  terrain  at  only  the  low  altitude 
end  of  the  link  is  considered.  Location  and  time  variabilities  are 
one-plecewist:  and  two-piecewise,  respectively,  normally  distributed  in 
both  programs. 

The  Longley-Rice  and  Johnson-Cierhart  programs  predict  the  hourly 

median  quantiles  L_  (d,q)  of  the  basic  transmission  L. (d)  as  a 
bo  b 

function  of  the  path  distance  d  and  the  cumulative  distribution 
function  (confidence  level)  q  =  prob  IL.  (d)  <  !  fe(j(d)  ] .  The  single 
variate  probability  q,  in  which  marginal  probabilities  of  time, 
location,  and  situation  (model  uncertainity )  are  combined,  is  the  only 
statistical  service  (designated  single-message  service)  available  in 
the  Johnson-Gie^hart  program.  In  the  Longley-Rice  program  other 
statistical  services  are  also  available  (see  section  3.0,  Table  5.) 


For  the  Johnson-Gierhart  program,  in  which  path  loss  Lfa(d)  is 
two-piecewise  normally  distributed  with  a  breakpoint  at  the  median 
value  ^^(d,  50*),  the  expected  value  (  L^(d)>  and  standard  deviation 

°L  (d)  are  given  by 
b 

<LK>=  Lk  (50*)  ♦  0.312  IL.  (90*)  ♦  L  (10*)  -  2LW  (50*)). 
b  bo  bo  bo  bo 

two-piecewise  normal  (2.5) 


{  5  ['-b,<90‘>  -  Lbo(50»»2*  5  1  1-b„l50‘>  -  L  bi10”12 

k  tLbo(90,)  ♦  LM(,0,)  -  2Lbbl50,,)2  !  'n  • 

two-piecewise  normal  (2.6) 


In  Eqs .  (2.5)  and  (2.6),  the  parameter  d  has  been  suppressed  in  order 
to  condense  the  notation. 
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For  the  Longley-Rlce  program,  In  which  path  loss  Lb(d)  for  the 

low-altitude  scenarios  Is  normally  distributed,  < Lh(d)>  and  o.  (d) 

D  Lb 

are  given  by 

<Lb(d)>*  Lb0(d,  501),  normal  (2.7) 

V">  =  vis  [Lbo<d'  90,1  -  Lb0td-  5°‘» 

D 

=  L  (d,  84.11)  -  L.  (d,  501),  normal  (2.8) 

DO  DO 

2qs.  (2.5)  and  (2.6)  reduce  to  Eqs .  (2.7)  and  (2.8)  respectively  for 

L.  (90!)  s  l.  (101). 
bo  bo 

Both  the  Longley-Rlce  and  Johnson-Gierhart  programs  utilize 
information  based  on  the  same  statistical-empirical  study  of  various 
terrain  profiles  and  propagation  measurements  conducted  primarily  In 
the  United  States.  In  that  study,  each  terrain  profile  was 
characterized  by  its  interdeclle  height  Ah  and  distance  to  the  radio 
horizon.  The  Interdeclle  height  is  the  difference  in  heights 
corresponding  ie  90!  and  10!  values  of  the  cumulative  distribution 
function  for  the  height  deviation  from  the  mean  surface  level.  When 
random  siting  Is  specified  for  an  antenna  above  a  surface  of  specified 
interdeclle  height,  the  programs  assign  a  median  distance  to  the  radio 
horizon  on  the  basis  of  information  derived  from  the 
statistical-empirical  study. 

The  statistical  parameters  Lb(d>  and  ^are  tabulated  in 
Table  1.  The  basic  transmission  loss  is  tabulated  for  Interdeclle 
terrain  heights  Ah  =  90  m,  225  m,  and  500  m  corresponding 
respectively  to  hills  (U.  S.  average  terrain),  mountains  (slightly 
shorter  than  those  in  the  Fulda  gap  of  Germany) ,  and  rugged  mountains 
(slightly  taller  than  those  in  Korea).  The  results  of  Table  1  ere 


for  average  ground  permittivity,  random  siting  of  the  antennas, 
atmospheric  refractivlty  at  the  earth*s  surface  equivalent  to  an 
effective  earth  radius  equal  to  4/3  earth's  geometric  radius, 
temporate  climate,  vertioal  polarization,  and  isotropic  antennas  in  a 
non-vegetati ve  open  area. 


The  expected  value  of  the  basic  transmission  loss  exceeds  the 

free-space  loss,  for  Ah  =  90  m,  225,  and  5C0  m,  by  1.1  dB,  0.7  dB, 

0.4  dB  respectively  for  the  high-altitude  scenario  no.  8a1  and  by 

32.5  dB,  3**.5  dB,  and  42.7  dB  respectively  for  the  low-altitude 

scenario  no.  9a.  The  standard  deviation  o,  ...  is  0  -  3  dB  for  the 

Lbld; 

high  altitude  scenarios  of  Table  1  and  10.9-  12.2  dB  for  the 
low-altitude  scenarios. 


Longley-Rice  predictions  of  basic  transmission  loss  for 
worldwide  environmental  conditions  are  tabulated  in  Table  2  for 
scenario  no.  9a  of  Table  1.  Variations  in  terrain  roughness  and 
surface  permittivity  have  an  appreciable  effect  on  the  expected  value 
of  transmission  loss.  Very  careful  siting  of  the  ground  antenna  can 
reduce  transmission  loss  by  approximately  6  rfB  over  that  for  random 
siting.  Variations  in  climate  and  atmospheric  refractivlty  have 
relatively  little  effect  on  transmission  loss  for  this  scenario.  In 
non-vegetated  areas,  vertically  polarized  waves  have  less  path  loss 
than  for  horizontal  polarization.  However,  in  vegetated  areas,  the 
reverse  may  be  true  (see  earlier  discussion  of  vegetative  loss  term). 


It  will  be  noted  that  transmission  loss  decreases  with 
increasing  terrain  roughness  for  terrain  roughness  less  than  or 
comparable  to  the  higher  antenna  height  above  the  ground  but 
increases  with  increasing  terrain  roughness  for  terrain  appreciably 
larger  than  the  higher  antenna  height.  The  reason  is  that 
obstructions  with  small  interdeolle  heights  do  not  appreciably  reduce 
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Table  2 


Basic  Transmission  Loss  for  Worldwide 
Environmental  Conditions 

Scenario  PtrtMCiri:  hj  •  2.5  a.  hj  -  61  m,  d  -  30  km.  f  •  50  KKz. 
non-vegetated  open  area 

*  Barallne  Envlronaantal  Parameters:  Ah  ■  90  a,  average  ground  permittivity, 
random  altlng,  K  "  4/3  earth  radlua,  tempo race  climate,  vertical  polarization. 


Environs 

Basic  Transmission  Lose 
Lb  (dB) 

<s> 

\ 

Terrain 

Roughness,  Ah  (m) 

■  ■ 

0  (Perfectly  Smooth) 

7.1 

5  (Water  or  Very  Smooth  Plalna) 

■rIM 

7.5 

30  (Slightly  Koilir.g  Plalna) 

HtlVB 

9.3 

60  (Rolling  Plains) 

10.3 

*90  (United  States  Average) 

128.5 

10.9 

225  (Mountains'* 

130.5 

11.5 

500  (Rugged  fountains)  1 

138.7 

12.1 

700  (Extremely 

Rugged  Mountains) 

145.5 

12.5 

Surface 

Permittivity 

c 

o(S/m) 

4 

.001 

(Poor  Ground) 

■  ' 

10.9 

*  15 

.005 

^Average  Ground) 

10.9 

25 

.02 

(Good  Ground) 

IIBTfJFHil 

10.9 

81 

5. 

(Sea  Watar) 

113.4 

10.9 

81 

.01 

(Fresh  Water) 

124.4 

10.9 

Siting  Criteria  For  Ground  Facility 

*  Random  Siting 

10.9 

Careful  Siting 

10.9 

Very  Careful  Siting 

■sni 

10.9 

Climate 

DM 

Equatorial 

10.9 

Contlnantal  Subtropical 

10.9 

Maritime  Subtropical 

10.9 

Deaert 

129.0 

10.9 

*  Contlnantal  Temporate 

128.5 

10.9 

Maritime  Temporate  Overland 

128.6 

10.8 

Maritime  Temporate  Oversea 

128.5 

10.9 

Atmospheric  Refractlvity 

mmm 

K  (Earth  Radius) 

Nb  (N-unita) 

1.23 

250 

1 

10.9 

* 

1.33 

301 

■  9 

10.9 

1.49 

350 

■mBi 

10.9 

1.77 

400 

127.2 

10.9 

Polarization 

*  Vertical 

128.5 

10.9 

Horizonte! 

133.2 

10.9 
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path  clearance  but  Instead  enhance  propagation  by  knife-edge 
diffraction.  For  obatruotions  with  large  lnterdeolle  heights,  the 
exponential  increase  of  the  path  loss  with  increasing  distanoe  Into 
the  shadow  region  of  the  obstructing  terrain  exceeds  any  reduction  in 
path  loss  obtained  by  the  knife  edges  of  the  obstructing  terrain. 

When  smooth  spherloal  earth  diffraction  Is  the  dominant  mode  of 
propagation,  surface  permittivity  and  particularly  ground 
conductivity  have  appreciable  effects  on  the  distribution  of  energy 
above  and  below  the  earth's  surfaoe.  For  example,  in  Table  2,  the 
transmission  loss  is  approximately  17  dB  less  for  propagation  paths 
over  sea  water  than  for  over  very  dry  (poor)  ground. 

For  non-vegetated  open  areas,  the  expeoted  value  and  standard 
deviation  of  the  excess  propagation  loss  A(d)  are  relatively 
frequency  Insensitive,  over  the  frequency  range  30  -  88  MHz,  when 
compared  to  the  variation  of  free-spaoe  loss  over  this  frequency 
range  (see  Table  3). 

A  comparison  of  theoretical  models  with  Longley-Rloe  predicted 
values  Is  given  In  Table  A  for  scenario  no.  9a  of  Table  1  and  an 
lnterdeclle  terrain  roughness  Ah  *  0.  The  theoretical  models  which 
are  considered  are  free  space,  plane  aarth  multipath,  and  spherloal 
earth  multipath.  The  seml-empirlcal  model  gives  an  expected  value 
which  exceeds  the  free-spaoe  loss  by  35.4  -  37.5  dB  over  the 
frequency  range  30  -  88  MHz.  The  plane-earth  multipath  model 
predicts  a  loss  whioh  Is  1  to  6  dB  larger  than  that  predicted  by  the 
Longley-Rlce  model.  The  spherical  earth  multipath  model  prediots  a 
loss  whioh  Is  appreciably  less  than  that  predloted  by  the 
Longley-Rloe  model  but  is  more  than  the  free-spaoe  loss.  The  close 
agreement,  between  the  results  for  the  multipath  plane  earth  model 
and  those  of  the  Longley-Rloe  model,  should  be  viewed  as  Just  a 
coincidence  because  the  multipath  plane  earth  model  is  only  an 


TABLE  3:  FREQUENCY  SCALING  0?  EXCESS  PROPAGATION  LOSS 
OVER  FREE- SPACE  LOSS 


ress 


idealized  limit  of  the  more  physical  spherloal  earth  multipath  model 
which  in  turn  is  inappropriate  for  radio  propagation  paths  that  clear 
the  radio  horizon  by  less  than  1/8  Fresnel  number.  The  Incremental 
path  loss  (expressed  as  a  numeric  rather  than  in  units  of  dB) 
increases  with  incremental  range  to  the  second,  fourth,  and  greater 
than  fifth  powers  for  the  free-space,  plane-earth  multipath,  and 
Table  2  Longley-Rlce  models  respectively. 

A  comparison,  of  the  Longley-Rlce  semi-empirical  model  with 

theoretical  models  and  with  empirical  data  in  a  multipath 

(17)-( 19) 

Interference  mode  of  propagation,  has  been  reported  .  It  was 

found  that  both  the  Longley-Rlce  semi-empirical  model  and  a 
statistical  model,  in  which  the  surface  height  is  assumed  to  be 
exponentially  distributed,  give  good  agreement  with  experimental  data 
for  coherent  scatter  in  the  forward-scatterd  direction  for  both 
terrain  and  sea  surfaces.  The  theoretical  model  for  the  particular 
mode  of  propagation  has  the  advantages  of  providing  a  theoretical 
basis  for  the  results  and  better  agreement  with  data  for  very  smooth 

surfaces  and  possible  very  rough  surfaces.  However,  for  the 

multipath  diffraction  transition  mode  of  propagation  near  the  radio 

horizon  over  an  irregular  terrain,  the  Longley-Rlce  semi -empirical 
model  appears  to  be  the  best  available  model  because  there  presently 
is  no  adequate  theoretical  model. 
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SECTION  3 


INPUT  PARAMETER  SPECIFICATION  FOR  LONGLEY-RICE,  VERSION  1.2.1, 

PREDICTION  PROGRAM 

The  Longley-Rice ,  version  1.2.1,  propagation  program  predicts 
long-term  (hourly)  median  radio  transmission  loss  over  irregular 
terrain.  The  output  of  the  program  is  basic  transmission  loss.  The 
program  combines  well-established  propagation  theory  with  empirical 
data  to  predict  propagation  lcsses.  The  prediction  program  is 
applicable  for  radio  frequencies  above  20  MHz.  For  frequencies  below 
40  MHz,  a  warning  is  automatically  printed  out,  regardless  of  path 
distance,  to  remind  the  user  that  the  sky  wave  may  be  significant  for 
sufficiently  long  paths.  The  program  may  be  used  either  with  terrain 
profiles  that  are  representative  of  median  terrain  characteristics 
for  a  given  area  (the  area  prediction  mode)  or  with  detailed  terrain 
profiles  for  actual  paths  (the  point-to-point  mode). 

The  empirical  data  base  is  for  wide  ranges  of  frequency,  antenna 
height  and  distance,  and  for  all  types  of  terrain  from  very  smooth 
plains  to  extremely  rugged  mountains.  The  data  base  includes  more 
than  500  long-term  recordings  at  fixed  locations  throughout  the  world 
in  the  frequency  range  40  MHz  to  10  GHz,  and  several  thousand  mobile 
recordings  in  the  United  States  at  frequencies  from  20  MHz  to  1  GHz. 
Much  of  the  empirical  data  base  is  in  the  VHF  frequency  band  30-100 
MHz. 


Preceding  Page  Blank 
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The  program  is  intended  for  use  within  the  following  ranges: 


Parameter 
frequency 
antenna  heights 
distance 


Range 

20  -  40,000  MHz 
0.5  -  3,000  m 
1  -  2,000  km 


surface  refractivity  250  -  400  N-units 

elevation  angle,  of  the  irregular  terrain 
radio  horizon  ray  above  the  horizontal, 

at  each  antenna  0-12  degrees 

relative  distance,  from  each  antenna  to  its 
terrain  horizon,  normalized  to  the 

corresponding  smooth-earth  distance  0.1  -  3.0 


The  elevation  angles  and  radio  horizon  distances  are  not  program  input 
parameters  but  are  computed  internally  by  the  program. 

Version  1.2.1  is  written  in  ANSI  Fortran  language  and  is 
therefore  compatible  with  any  large  scale  computer. 

The  input  parameter  specifications  for  version  1.2.1  are  given  in 
Table  5  which  includes  numerical  ranges  of  parameters  and  the 
numerical  values  for  which  a  warning  is  automatically  printed  out.  An 
asterisk  denotes  the  numerical  value  that  will  be  assumed  for  a 
parameter  if  the  user  does  not  specify  a  particular  value. 

Version  1.2.1  offers  two  program  modes  whose  selection  depends 
upon  how  the  user  wishes  to  specify  the  terrain  surface  profile:  the 

area  prediction  mode,  designated  "QKAREA";  and  the  point-to-point 
mode,  designated  "QKPFL".  The  area  prediction  mode  is  characterized 
by  specifying  the  interdecile  height  ah,  the  antenna  siting  criteria, 
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INPUT  PARAMETER  SPECIFICATION,  LONGLEY-RICE  PREDICTION  PROGRAM,  VERSION 


fAMCTd 


TABLE  5:  INPUT  PARAMETER  SPECIFICATION,  LONGLEY-RICE  PREDICTION  PROGRAM,  VERSION  1.2.1  (con 


TABLE  5:  INPUT  PARAMETER  SPECIFICATION.  LONOLEY-RICE  PREDICTION  PROGRAM,  VERSION  1.2.1  (cont. 


and  the  great  circle  distances  for  which  the  basic  transmission  loss 

L. (d)  is  printed  out.  The  point-to-point  mode  is  characterized  by 
b 

specifying  the  elevation  matrix  of  the  terrain  profile.  In  the  area 
prediction  mode,  the  lnterdeclle  terrain  height  and  antenna  siting 
criteria  determine  the  expected  values  of  the  antenna  effective 
heights,  the  elevation  angles,  and  the  terrain  horizon  distances  by 
means  of  the  stored  empirical  data  base.  The  point-to-point  mode 
computes  the  elevation  angles  and  distance  to  the  radio  horizon  for 
each  antenna  and  the  antenna  effective  height  by  considering  whether 
or  not  the  antenna  is  near  or  on  a  hill.  Except  for  these 
differences,  the  QKAREA  and  QKPFL  programs  are  identical. 

The  in out  card  set-ups  lor  programs  QKAREA  and  QKPFL  are  given  in 
Tables  A-l  and  A-2  of  Appendix  A,  respectively.  The  parameter  fields 
for  programs  QKAREA,  QKPFL  and  program  control  cards  are  given  in 
Tables  A-3,  A-*»  and  A-S  of  Appendix  A,  respectively. 

A  card  deck  consists  of  Job  control  cards,  input  parameter  card 
types,  and  program  control  card  types. 

There  are  four  different  job  control  cards.  They  preceed  the 
rest  of  the  deck  and  all  are  required  in  the  sequence  shown  in  the 
card  set-up  tables. 

Input  parameter  card  types  1  -  7  follow  the  Job  control  cards. 

For  QKAREA,  there  are  as  many  as  six  card  types  for  each  "execute" 
operation.  For  QKPFL,  there  are  as  many  as  seven  card  types  for  each 
"execute"  operation.  Any  parameters  not  specified  (left  blank) 
default  to  values  noted  in  the  input  parameter  specification  table. 

Any  card  type  can  be  omitted.  If  a  card  type  is  omitted,  the  program 
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defaults  to  the  values  given  In  the  input  parameter  specification 
table.  If  a  card  type  Is  used,  all  cards  In  that  card  type  must  be 
specified  In  the  sequence  shown  in  the  card  set-up  tables.  Card  types 
1  -  7  can  be  Inserted  In  any  order. 

There  are  three  program  control  card  types:  0  (stop),  8 
(execute)  and  9  (reset).  The  functions  of  these  cards  are: 

stop  -  causes  the  program  to  execute,  produce  a  printout  and 
then  terminates  the  Job. 

execute  -  causes  the  program  to  execute  and  produce  a  printout 
for  a  given  set  of  Input  parameter  card  types. 

reset  -  causes  all  parameters  to  be  reset  to  their  default 
values  until  additional  input  parameter  card  types  are 
inserted . 

The  program  control  cards  are  placed  after  the  four  Job  control  cards 
at  approximate  places  within  the  card  deck  where  their  functions  are 
required.  Card  type  0  is  required  for  final  execution  and  should  be 
the  last  card  in  the  deck.  Card  type  8  is  placed  in  the  card  deck 
before  each  set  of  input  parameter  cards  types.  Card  type  9  is  used 
only  if  there  is  more  than  one  set  of  input  parameter  card  types 
required  for  additional  runs  and  is  placed  after  the  execute  card  but 
before  each  additional  set  of  input  parameter  cards. 

Sample  outputs  for  the  longley-Rice  QKAREA  and  QKPFL  programs  are 
given  in  Tables  6  and  7  respectively. 
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TABLE  6 


Sample  Output,  QKAREA 


PROCRAM  QKAREA 


PREQUENCT 

80. 

,  MHZ 

ANTENNA  HEIGHTS 

2.8 

78.0  M 

EFFECTIVE  HEIGHTS 

2.9 

76.0  M 

TERRAIN,  DELTA  H 

700.  M 

(SIT INQ*0, 0 ) 


POLO >  EPS* 1 S. .  SOM*  O.OOS  S/M 
CLlM.fc,  NO ■ 30 1 • <  NS*S01..  K*  1.333 


SINGLE-MESSAGE  SERVICE 

ESTIMATED  QUANTILES  OP  SASIC  TRANSMISSION  1058(09) 


OIST 

FREE 

HI  Til  CONFIDENCE 

50.0 

70.0 

KM 

SPACE 

5.0 

16.9 

5.0 

80.4 

06.4 

94.8 

107.2 

113.6 

10-0 

86.4 

96.6 

105.0 

117.3 

123.8 

lS.0 

90.0 

104.3 

<12.3 

125.0 

131 .4 

20.0 

92.8 

111.0 

119.3 

131 .6 

138.0 

25.0 

94.4 

117.1 

126.4 

137.7 

144.1 

30.0 

96.0 

122.8 

131.3 

143.4 

149.8 

35.0 

97.3 

128.3 

136.6 

148.9 

155.2 

49-0 

96. 5 

133.6 

141.9 

154.2 

160.5 

45.0 

99.5 

137.0 

149.4 

157.6 

163.9 

SO  .0 

100.4 

139.0 

147.4 

159.7 

166.0 

55.0 

101 . 3 

140.9 

149.3 

181 .7 

168.0 

60.0 

102.0 

142.6 

181 .2 

163.6 

169.6 

70.0 

103.3 

145.8 

184.6 

187.1 

173.4 

80.0 

104.6 

148.6 

157.7 

170.8 

176.9 

90.0 

105.6 

161  .2 

160.6 

173.7 

180.1 

100.0 

109.4 

183.9 

163  3 

176.8 

1C3.3 

110.0 

107.3 

153.6 

169.9 

179.8 

ice.  3 

120.0 

108.0 

198.0 

168.4 

163.9 

189.3 

84.1 

119.  S 

139.7 

137.4 
144.0 
150.0 

155.7 

191.1 

166.3 

169.6 
171  .8 

173.6 

178.7 

179.3 

103.7 

186.1 

189.3 

193.4 

198.4 


90.0 

133.3 

133.3 
141.0 

147.8 

183.6 

159.3 

164.6 

169.9 

173.3 

175.4 

177.3 

179.3 

183.9 

166.3 

169.7 
193.0 

196.1 

199.2 


99.0 

137.9 

136.1 

145.7 

193.3 

168.3 

103.9 

169.3 

174.8 
178.0 

180.1 
163. 0 

193.9 
187.6 

191 .1 

194.8 

197.9 

201.1 

204.3 
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TABLE  7:  Sample  Output,  QKPPL 


g lift  TEST  1,  PATH  2200  (MEASURED  MEDIAN  LB-132, 2  DB) 
CRYSTAL  PALACE  TO  MURSLEY,  ENGLAND 


DISTANCE 
FREQUENCY 
ANTENNA  HEIGHTS 
EFFECTIVE  HEIGHTS 
TERRAIN,  DELTA  H 


77.8  *KM 
41.5  MHZ 
143.9  8.5  M 

240.5  18.4  M 

89.  M 


POL-O ,  EPS-15.,  8GM-  .005  S/M 
CLIM-5,  N8-314. ,  K*  1.368 
PROFILE-  HP-  156,  XI-  .499  KM 


A  DOUBLE-HORIZON  PATH 
DIFFRACTION  IS  THE  DOMINANT  MODE 


ESTIMATED  QOANTZLE8  OF  BASIC  TRANSMISSION  LOSS  (DB) 
FREE  SPACE  VALUE-  102.6  DB 


RELIA¬ 

WITH  CONFIDENCE 

BILITY 

50.0 

90.0 

10.0 

1.0 

128.6 

137.6 

119.6 

10.0 

132.2 

140.8 

123.5 

50.0 

135.8 

144.3 

127.2 

90.0 

138.0 

146.5 

129.4 

99.0 

139.7 

148.4 

131.0 
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SECTION  U 


INPUT  PARAMETER  SPECIFICATION  FOR  JOHNSON-GIERHART , 
AIR-TO-AIR  ( ATOA ) ,  PREDICTION  PROGRAM 


The  Johnson-Gierhart  ATOA  prediction  program  is  similar  in  many 
respects  to  the  Longley-Rice  prediction  program.  The  ATOA  program 
predicts  radio  transmission  loss  over  irregular  terrain.  The  output 
of  the  MITRE  versions  of  the  program  is  basic  transmission  loss.  The 
program  combines  well-established  propagation  theory  with  empirical 
data  to  predict  propagation  losses  which  have  been  tested  against  a 
large  number  of  propagation  measurements.  It  is  applicable  to  radio 
frequencies  100  MHz  to  20  GHz  but  can  also  be  used  at  frequencies  as 
low  as  20  MHz  provided  that  the  propagation  paths  are  sufficiently 
short  so  thut  the  sky  wave  is  insignificant.  The  program  is  used  with 
terrain  profiles  that  are  representative  of  median  terrain 
characteristics  for  a  given  area  (similar  to  the  area  prediction  mode 
of  version  1.2.1  of  the  Longley-Rice  program)  but  not  with  detailed 
terrain  profiles  for  actual  paths. 

The  Johnson-Gierhart  ATOA  prediction  program  differs  from  the 
Longley-Rloe  program  principally  in  the  following  ways: 

a)  Only  single-horizon  diffraction  rather  than  double-horizon 
diffraction  is  considered. 

b)  A  standard  exponential  atmosphere  rather  than  a  uniform 
gradient  atmosphere  is  assumed  (the  index  of  refraction 
decreases  exponentially  rather  than  linearly  with  increasing 
height) . 
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c)  Plotting  routines  and  various  output  options  are  available 
in  the  original  CDC  computer  versions  of  the  ATOA  program 
but  are  presently  not  available  in  the  ATOA  MITRE  versions. 

d)  Basic  transmission  loss  options  exist  for  specifying  whether 
surface  reflection  multipath  and  tropospheric  multipath 
contribute  to  Instantaneous  levels  exceeded  or  hourly  median 
levels  exceeded  and  whether  It  contributes  to  variability 
(standard  deviation)  or  median  level  (501  confidence  level). 

Thj  empirical  data  base  comprises  200  single-horizon  paths  from 
Reference  [16].  Double-horizon  paths  are  not  included  in  the 
empirical  data  base. 

The  program  is  intended  for  use  within  the  following  ranges: 


Parameter 

frequency 

lower  antenna  height 
higher  anenna  height 


surfaoe  refractlvity 

elevation  angle,  of  the  Irregular  terrain 
radio  horizon  ray  above  the  horizontal, 
at  the  lower  antenna  only 

distance,  from  the  lower  antenna  to  its 
terrain  horizon,  relative  to  the 
corresponding  smooth-earth  distance 


Range 

100  -  20,000  MHz 

0.5  —  3*000  m 

radio  horizon 
height  of  lower 
antenna 

200  -  «00  H-units 


0-12  degrees 


0.1  -  3.0 
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The  input  parameter  specifications  for  the  ATO A  program  are 
given  in  Table  8.  An  asterisk  denotes  the  numerical  value  that  will 
be  assumed  for  a  parameter  by  the  progranmer  if  the  user  does  not 
specify  a  particular  value.  Unlike  the  Longley-Rioe  version  1.2.1 
program,  the  asterisk  values  are  not  program  automatic  default  values 
but  must  be  specified  by  the  programner .  The  elevation  angle  of  the 
horizon  at  the  facility  (lower  antenna)  and  •the  distance  from  the 
lower  antenna  to  its  radio  horizon  are  not  program  input  parameters 
but  are  calculated  Internally  by  the  program. 

The  ATOA  program  which  was  originally  written  for  a  CDC  6*1  bit 
word  computer  has  been  converted  by  MITRE  for  use  on  its  IBM  370  32 
bit  word  computer.  The  converted  ATOA  source  module  has  been 
compiled  on  the  IBM  03  FORTRAN  IV  H  Extended  Compiler  using  the 
AUTOBBL  (DBLPAD)  and  OPTIMIZE  =  0  options. 

The  converted  program  contains  only  the  output  option  of  basic 
transmission  loss  and  has  neither  the  plotting  capability  nor  the 
other  output  options  of  the  original  program.  The  converted  program 
gives  results  which  are  identical  to  those  of  the  original  test  run 
dated  04/25/79. 

MITRE  has  three  versions  of  ATOA.  Version  1  produces  numerical 
output  which  is  identical  to  the  original  except  for  a  few  minor 
differences  in  format.  Version  2  allows  the  option  of  speolfylng 
confidence  levels,  distance  increment  values,  and  a  condensed 
parameter  heading  on  print  out.  Version  3  has  the  same  options  as 
version  2  but  also  prints  out  the  expected  value  and  standard 
deviation  of  basic  transmission  loss  at  eaoh  distance. 


psihari  parameters,  specification  required 


Value 


TABLE  8:  INPUT  PARAMETER  SPECIFICATION,  JOHNSON-CIERHART  ATOA  PREDICTION  PROGRAM  (cont. 


The  Input  card  set-up  for  these  three  MITRE  versions  of  ATOA  is 
given  in  Table  B-1  of  Appendis  B.  The  parameter  fields  for  these 
versions  are  given  In  Table  B-2.  Important  files  related  to  MITRE 
versions  of  ATOA  are  listed  and  described  in  Table  B-3. 

A  sample  output  of  the  ATOA  program,  MITRE  version  3,  Is  given 
In  Table  9.  The  sample  output  corresponds  to  the  test  run  dated 
04/25/79  of  the  original  ATOA  program. 
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TABLE  9:  Sample  Output,  ATOA  MITRE  Version 

Mi  ATOA  (MI TAC  VERSION  3)  PACE  1  05/06/82  09.56-0* 
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S  %  ®oo-nflf  t«e#nfuiN#o--fn,Hf  ttf  <»«i»ib«md(|ioo 

—  -  ftncifinnpirinfiiinnricicinwfln 

5  s . - . - . ------------------- 

5  «« 

4  v  ;nnnn*invi5«o«»in»#o»--n«t*f<iinlfl*f#»nnnninni<»n,i»«ift*(Oh®» 

■  oooooooodO'-'o'-rinnNftNnnrinnnPinNnNM'inrinnnnnnnrtnnn 

-  | . . . . . .  ------ . - 

u. 

O  g  •-  - 

IA  O  . . . . **. . . . . 

v  o  a  m  4  a  a  n*  a  *  a  p*  a  r*  a  -  a  n  <*  «  a  a  e»  r*  f  r-  •  r-  »■*  a  a  a  •  *  a  n  a  a  n  n  4  m  m  ir  a  a  ^  a  » 
v  Oftoooooooo----NNnnnftnnnnNNfictnNnnftftnfinflrinnnnAnnMn 

j  S  . 

3  S  o  n  a  e  a  n  ®  a  n  n  a  e  a  <a  a  a  ia  a  p  a®  a  ^  a  a  -  *».  o  a  a  ®  m  n  ®  n  a  m  •  ir  m  ®  «*  a  o  At  n 

2  i  ^^§55^^5555 

► 

s  j  **xs;t***!!55$?!HHn?r!?iH*!!?-555555555555555555 

2  I«XSi«S*88St5»5:25:SS!52t5;;SS885S55S5S5S3SaR25* 

! 

.  §ssr?:ss£st;s;*3c*£si*8jj;i^5ts£8ssxssssss2sa*ss 
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u» 


TABLE  9:  Sample  Output,  ATOA  MITRE  Version  3  (CONT. ) 


>ooooooooo«ooo 


a<*ooooooooo  xSSSSsISSiixSxSS 


it*  i 

NNld 

••-a  * 

m  ? 

»  w 

Si*** 

^  —  to*  *  a 
•--auiO 

<  «  ►  ►  Ml 

22--S 

llitfwe 

1Z%*$ 

OOk-*- 

«*  *  ~ 
u  o  ¥  * 
?  a  w  £ 
2*  8  • 


s;ssssRR8saa5?t*R;?;js5sss52!esssg;s3ss*8S 


•  al  *i  o  o  o  -'  ti«i «  V  ir  n  o ft  V  *  h  «  » r*e-  o  e  «V  «  ^ «  f«  *  -'  *  m •  fli n  » 
—  -  •>rtdr<tfMr«c«drif^dddnnnnndm’«tft/)2^J»^4vrft^«2w^isr>d 


cca«aa999>»( 


id«0«d»o««r>^0>-«v 


•  nm  «  r>  «  r«  « 


n«h«^f  #a»»e«i»«A«rir*r»MrnaB*caMqhno»«oo 

nnnrtfih«««««««ui^i«tf2«hhN«aSo6-NAftnNnnnnnn 

^^^...rrfiveeeeffreveerfr^fMNmMnnNRflftnnnn 


III  : 


§"■  > 
w  o 

as- 

SO  S  tfl 

«-  4  w 

:*►  5 

si  .ss 

•  **l 

is  r. 

i?S..5 

55"S8 

liifi 

*s<rfd»u 
•*■  4CM 

SS23S 


liiili 


nnnflnnn<««<««^)AiA«ii 


'  O  O  -  - - -  •■  -  -  -  a*  I 


NNnnnn8nn'in<tf««o«iit#«e4«NM«««{»»A«iftSaoo 

sssaascRssRjsssssssssj^csssgsssstsssgss 

**^^*^»**e»er»er*a«a-ee«eeeeaaaaW^*-f|N(|»l 


r55rr5=5?5r5;55r5::???:r????H?????^E:?? 

cRfiRRassaaR**? 

55??!55555??5!?5^?5?55555!*555?!??^5? 

m:s!t;sssiui:iit»«s}tmms!sssssss 
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APPENDIX  A 


CARD  SET-UP  FOR  LONGLEY-RICE  QKAREA 
AND  OCPFL  PROGRAMS 


TABLE  A-2.  Card  Set-Up,  QKPFL 


.(Every  five  coluans  contains  one  elevation  P.  •  05  k  <n, 
continuing  on  to  as  many  cards  as  necessary. 7 


Up,  QKPFL 


53 


TABLE  A-3.  PARAMETER  FIELDS,  QKAREA  (cont.) 


TABLE  A-3.  PARAMETER  FIELDS,  QKAREA  (cont.) 


Card  Type  2 


Columns 

Algebraic 

Symbol 

Compucer 

Symbol 

Integer  (I) 
or  Floating 
Point  (F) 

Description 

1 

None 

None 

— 

Left  blank 

2 

None 

None 

I 

The  integer  2 

3-10 

Noae 

None 

— 

Left  blank 

11-20 

d 

0 

Vt 

F 

Great  circle  distance  (ka)  between 
terminals  (Initial  value). 

21-30 

di 

Dl 

F 

Great  circle  distance  (ka)  between 
terminals  (final  value,  in 
steps  of  dgl). 

31-40 

d.i 

DS1 

F 

Stepping  distance  1  (ka) . 

41-50 

d2 

D2 

F 

Great  circle  distance  (ka)  between 
terminals  (final  value.  It.  -teps 
of  d,2). 

51-60 

“.2 

DS2 

F 

Stepping  distance  2  (ka) . 

61-80 


None 


None 


Left  blank 


TABLE  A-3.  PARAMETER  FIELDS,  QKAREA  (cont.) 


C«rd  Type  3 


Symbol 

r.v*mpuc«r 

Symbol 

Integer  (l> 
or  '’looting 
Point  (F) 

je<cr  Ipt ion 

i 

None 

Non* 

Left  blank 

> 

None 

Non* 

t 

Th«  Integer  3 

i 

V 

XDVAJt 

i 

Mode  of  variability: 

(0)  -  single-message  service 

(1)  -  individual  service 

(2)  -  mobile  service 

(3)  -  broadcast  service 

i:-:o 

qT  or  q. 

QT  or  QR 

F 

Tine  or  location  reliability  (1) 

11-30 

^L 

QL 

f 

Location  reliability  (1) 

31-80 

|  Nor.e 

Non* 

1 

Left  blank 

TABLE  A-3.  PARAMETER  FIELDS,  QKAREA  (cont.) 

Card  Type  4 


Algebraic 
Symbol 


Integer  (I) 
Computer  or  floating 

Symhol  Point  (F) 


Descript  Ion 


Left  black 

The  integer 

4 

Laft  blank 

Conf ldenca 

level 

1 

(1) 

Confidence 

level 

2 

(X) 

Confidence 

level 

3 

(X) 

Conflcence 

level 

4 

(X) 

Confidence 

level 

5 

(X) 

Confidence 

level 

6 

(X) 

Confidence 

level 

7 

(X) 

SUc. 


TABLE  A-3.  PARAMETER  FIELDS,  QKAREA  (cont.) 


Card  Type  5 


Column* 

Algebraic 

Symbol 

Computer 

Symbol 

Integer  (1) 
or  Floating 
Point  <F) 

Dcscript Ion 

a'  — 1 

i 

Nona 

-- 

Left  blank 

2 

Nona 

I 

The  integer  S 

3 

C 

I 

Climate  code: 

(1)  -  equatorial 

(2)  -  continental 

(3)  -  maritime  subtropical 

(4)  -  desert 

(5)  -  continental  temperate 

<6)  -  maritime  temperate  overland 
(7)  -  maritime  temperate  oversea 

4-10 

None 

.. 

Lefc  blank 

11-20 

Ah 

DELTAH 

F 

Terrain  Irregularity  (m) 

21-30 

or 

N 

o 

EN0 

F 

Minimus  monthly  mean  of  atmospheric 
refractlvlty  at  sea  level. 

21-30* 

N 

8 

ENS 

F 

Atmospheric  refractlvlty  at  average 
elevation  of  ground  surface. 

31-40 

z 

9 

ZSYS 

1 

F 

1 

Average  elevation  of  ground  surface 
above  naan  sea  level  (a) . 

41-50 

c 

EPS 

F 

Surface  dielectric  constant. 

51-60 

0 

SCH 

F 

Surface  conductivity  (S/n) . 

61-80 

None 

1 

Nona 

1 

Left  blank 

*N  can  b 
8 

J  specif!  as  In  col 

Istead  of  N  . 

0 

60 


4-.J5.-WW~* 


TABLE  A-3.  PARAMETER  FIELDS,  QKAREA  (cont.) 


Card  Type  6 


Co  1  unins 

Algebraic 

Symbol 

Computer 

Symbol 

Integer  (l) 
or  Floating 
Point  (F) 

Descript  ion 

i 

Non* 

None 

— 

Left  blenk 

2 

None 

None 

I 

The  Integer  6 

3 

None 

None 

— 

Left  blenk 

U 

P 

POL 

I 

Antenna  polarization 

(0)  -  horizontal 
(1)  -  vertical 

5 

S1 

KST1 

X 

Siting  criteria  for  terminal  1 

(0)  -  random  siting 

(1)  -  careful  siting 

(2)  -  very  careful  siting 

6 

S2 

KST2 

I 

Siting  criteria  for  terminal  2 
(same  options  as  above). 

7-10 

None 

None 

— 

Blank 

11-20 

f 

FMHZ 

F 

Frequency  (MHz) 

21-30 

h  i 

gl 

HG1 

? 

Antenna  height  (■)  above  ground 
at  terminal  1. 

31-40 

h  , 
g^ 

HC2 

F 

Antenna  height  (m)  above  ground 
at  terminal  2. 

41 -BO 

None 

Non* 

Left  blank 

61 
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TABLE  A-4 .  PARAMETER  FIELDS,  QKPFL 


Card  Type  1 
(CarJ  1.1) 


1 

Integer  (I) 

Algebraic 

Computer 

or  Floating 

Columns 

Symbol 

Symbol 

Point  (F) 

Descript  ton 

TABLE  A-4 .  PARAMETER  FIELDS,  QKPFL  (cont.) 


1 


1 

\ 


Card  Type  1 
(Card  1.2) 


Integer  (I) 

Algebraic 

Computer 

or  Floating 

Columns 

Symbol 

Symbol 

Point  <P) 

Oeacrlptlon 

1-60 

None 

Nona 

Alpha/Niaerlc 

Output  title 

61-80 


Nona 


Nona 


Laft  blank 


TABLE  A~4.  PARAMETER  FIELDS,  QKPFL  (cont.) 


Card  Type  2 
(Card  2.2) 


(ulumn* 


Algebraic 

Symbol 


Computer 

Symbol 


Integer  (t) 
or  Floating 
Point  (F) 


Description 


66 


QKPFL  (cone.) 


Descript  ion 


This  coltan  contains  tbs  lotsgar  1 
If  this  card  la  tba  last  (or  only) 
profile  card.  Otherwise,  It  Is  left 
blank. 

These  colusna  contain  the  number  of 
profile  points  specified  on  this 
card.  (1  or  2  digits,  right 
Justified) 

Left  blank 


Profile  Matrix  (meters). 

(A  decimal  point  may  be  substituted 
for  one  of  the  Integers  within  the 
flvs  column  profile  fields,  if 
a  decimal  point  Is  not  specified,  s 
decimal  point  will  be  assrned 
after  the  last  lntagsr  in  thst 
field.) 


TABLE  A-4.  PARAMETER  FIELDS,  QKPFL  (cont.) 


Card  Type  3 


Co  luf  u  »3 

Algebraic 

Symbol 

Computer 

Syabol 

Integer  (I) 
or  Floating 
Point  (F) 

Description 

i 

None 

None 

Left  blank 

2 

None 

None 

I 

The  Integer  3 

3-10 

None 

None 

— 

Left  blank 

11-20 

qTl 

QTl 

t  \ 

21-30 

qT2 

QT2 

F 

31-40 

qT3 

QT3 

F 

41-50 

qT4 

QT4 

F 

Relleblllty  eetrlx  (X) 

51-60 

qT5 

QT5 

F  | 

, 

61-70 

qT6 

QT6 

T 

71-80 

i  <1T7 

QT7 

7  7 

I 


I 


r 


TABLE  A- A. 

PARAMETER  FIELDS,  QKPFL  (cont,) 

Catd  Type 

4 

I  Algebraic 

j  Computer 

I  Integer  (I)  1 

or  Floating  ! 

Columns 

|  Symbol 

1  Symbol 

|  Point  (F)  ]  Description 

Left  bunk 
The  integer  4 
Left  bUnk 


Confidence  level  matrix  (I) 


TABLE  A- 4 


PARAMETER  FIELDS,  QKPFL  (cont.) 


Card  Type  7 


Columns 

Algebraic 

Symbol 

Computer 

Symbol 

Integer  (1) 
or  Fleeting 
Point  (F) 

Description 

i 

None 

None 

Lefc  blank 

2 

None 

None 

I 

The  integer  7 

3 

P 

P01. 

I 

Antenna  polarliation: 

(0)  -  horlxoncal 
(1)  -  vertical 

u 

C 

aiM 

I 

Climate  code 

(1)  -  equatorial 

(2)  -  continental  aubtroplcal 

(3)  -  ear it lee  subtropical 

(4)  -  desert 

(5)  -  continental  temperate 

(6)  -  maritime  temperate  overland 

(7)  -  maritime  temperate  oversea 

5-10 

None 

None 

-- 

Left  blank 

11-20 

f 

FMHZ 

F 

Frequency  (MHt) 

21-30 

h  i 

8> 

HG1 

F 

Haight  above  ground  (m)  at  terminal  l. 

31-40 

h  •> 

g2 

MG2 

F 

Height  above  ground  (m)  at  terminal  2. 

o 

1 

w* 

o 

No 

EN0 

F 

Minimum  monthly  m-an  of  atmospheric 
refractlvlty  at  sea  level. 

51-60 

\ 

ENS 

F 

Atmospheric  refractlvlty  at  average 
elevation  of  ground  surface. 

61-70 

c 

EPS 

F 

Dielectric  Constant 

71-80 

0 

sen 

F 

Surface  conductivity  (S/m) 

69 


TABLE  A-5.  PARAMETER  FIELDS,  PROGRAM 
QKAREA  AMD  AKPFL 

Card  Type  8 


Column* 


Algebraic 

Symbol 


Computer 

Symbol 


Integer  (I) 
or  Floating 
Point  (F) 


CARDS 


■  nmMtf  ■ 


1 

r* 

% 


4 

■y 


TABLE  A-5.  PARAMETER  FIELDS,  PROGRAM  CONTROL  CARDS, 
Q KARRA  AND  QKPFL  (cont.) 


I 

L 


Card  Typa  0 


Column 8 

1 

2 


Algebraic 

Symbol 


None 

None 


Coapucer 

Symbol 


None 

None 


Integer  (1) 
or  Floating 
Point  (F) 


Description 


Left  blank 
the  Integer  0 


1 


3-80 


None 


None 


Left  blank 


TABLE  A- 5 ,  PARAMETER  FIELDS,  PROGRAM  CONTROL  CARDS, 
QKAREA  AMD  QKPFL  (cont . ) 


Columns 

Algebraic 

Symbol 

Computer 

Symbol 

Integer  (1) 
or  Floating 
Point  (F) 

Description 

1 

None 

Nona 

Left  blank 

2 

None 

None 

t 

The  lntagar  9 

3-80 


None 


None 


Left  blank 


APPENDIX  B 


CARD  SET-UP  FOR  JOHNSON -G I ERH ART  ATOA  PROGRAM , 
MITRE  VERSIONS  1,2  AND  3 


73 


TABLE  B-2 .  PARAKETER  FIELDS,  ATOA 


MITRE  VERSIONS  1,  2  &  3 
CARO  1 


Card  | 

Columns' 

Symbol 

Description 

- - l_ 

i 

l'2  1 

IK 

Code  for  units  to  be  used  with  input: 

! 

i 

(1)  las  and  meters 

(2)  ft  and  st  ml 

(3)  ft  and  nml 

(4)  IK  ■  0  terminates  a  run 

3-4  | 

10 

Code  for  type  of  output: 

! 

(1)  Power  available 

(2)  Power  Density 

(3)  Transmission  Loss 

5-6 

j 

1J 

Code  for  aircraft  altitude  input: 

1 

>0  the  units  will  be  considered  the  same  as  distance, 

1 .e.  km,  a  ml  or  n  ml . 

7-8 

I  LB 

Code  for  loblng  options: 

1 

(0)  No  loblng,  (2)  lobing. 

9-10 

KK 

Code  for  time  availability  options: 

(1)  hourly  median  levels,  (2)  instantaneous  levels 

11-13 

1A 

Number  of  characters  and  spaces  in  label. 

May  be  up  to  32.  (Not  used  in  MITRE  versions) 

14-45 

TT 

Label 

77 


TABLE  B-2.  PARAMETER  FIELDS,  ATOA  (cont.) 


:;im  ver*:  >ns  i .  2  a  3 

CAKD  2 


Card  j 


Columns! 

Symbol 

Description 

-  1 

1 

:-s  i 

i 

DMIN 

1  Abscissa  value  for  left-hand  limit  of  graph  (n  mi  or  deg.), 
(used  by  ICPH) 

6-10  j 

1 

DKAa 

Abscissa  value  fnr  right-hand  limit  of  graph  (n  ml  or  deg.). 

|  (used  by  ICPH) 

1  i  —  1 4 

1 

t 

1 

xc 

Abscises  increment  for  graph  grid  lines  (n  ml  or  deg.). 

(Not  used) 

15-19  j 

P.-liN  j 

1 

1 

[ordinate  value  for  bottom  limit  of  graph: 

(dB-W/sq.  ml  for  Power  Density,  dBU  for  Power 
|  Available  and  dB  for  Transmission  Loss)  (Not  used) 

I 

Must  be  negative  for  Transmission  Loss. 

20-24  j 

PMAX  | 

I  Or  Unate  value  for  top  limit  of  graph 

i 

i 

i 

Must  be  negative  for  Transmission  Loss. 

25-28 

YC 

Orllnate  Increment  for  graph  grid  lines  (Not  used) 

Must  be  positive  for  Transmission  Loss. 

2  7-30  | 

JC 

I 

If  the  output  Is  to  be  plotted  against  deg  JC  >  0. 

31-32 

ICPH  ' 

Code:  <1  Plotting:  (t)  No  Plotting;  (>1)  Will  Interpolate 

and  get  values  for  distance  In  DMAX,  columns 

51-55.  Also  no  Plotting. 

(CODE:  (0)  not  available). 
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TABLE  B-2,  PARAMETER  FIELDS,  ATOA  (cont.) 


illTRl  VERSIONS  1.243 
CARD  3 


Card  | 
Columns! 


Deacrlptlon 


Height  of  facility  or  lower  antenna  above  mil 
Code  for  facility  antenna  pattern: 

(1)  Isotropic 

(2)  DME 

(31  TACAN  (RTA-2) 

(4)  4-loop  array  (cosine  vertical  pattern) 

(5)  8-loop  array  (cosine  vertical  pattern) 

(6)  1  or  It  (cosine  vertical  pattern) 

(7)  JTAC  with  tilted  antenna 

(8-21)  Special  antennas 


10-11 

!  -,T  | 

Code  for 

antenna:  (0)  directive  (»0)  tracking 

12-13 

1  ipl  i 

Code  for 

polarization  of  facility  antenna: 

j  | 

i 

(1) 

Horizontal 

;  i 

(2) 

Vertical 

:  | 

(3) 

Clrcul ar 

1 

Also 

used  for  reflection  coefficient  and  ground  constants 

14-18 

i  TIT 

Tilt  of 

the  facility  antenna  naln  beam  In  deg. 

|  Not 

used  for  Patterns  1-6. 

19-23 

'  H1.PBW  j 

jHalf  of 

the  lialf-pouer-beaa  width  of  the  facility  antenna 

1 

!'  Nor 

1 

used  for  patterns  1-e. 

24-28 

SUR 

I 

Elevation  of  facility  site  surface  above  msl 

29-30 

1 

IZ 

Rainfall 

Zones  (0)  no  consideration  (1-6)  see  Samson's  maps. 

1  (?) 

Adds  0.3  dB  times  storm  size  to  attenuation. 

31-33 

;  STS 

\  Size  of 

storm:  3,  10,  or  20  km. 

34-33 

|  KD 

1 

iCode  for 

terrain  type  options: 

(1) 

smooth  earth 

(2) 

Irregular  terrain. 

36-37 

K£ 

Code  for 

horizon  options: 

(0) 

none  specified 

(1) 

angle  specified  by  1  DC ,  IMN,  and  SEC 

(2) 

height  specified  by  HH01 

(3) 

both  the  any, ip  and  the  elevation  are  specified. 

38-4  1 

UHS1 

Tet  rain 

parameter  6h  (ft)  from  cable. 
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TABLE  B-2.  PARAMETER  FIELDS,  ATOA  (cont.) 


MITRE  VERSIONS  1.  2  i  3 
CARD  4 


Card  i 

Coluuny^  Symbol 


28-33  1 


Description 


DHOI  Distance  to  facility  radio  horizon  (n  ml). 

NOTE:  Zero  or  negative  values  will  result  In 
calculation  of  this  parameter  from  othera  (fig.  14). 

HHOI  Elevation  of  facility  radio  horizon  above  ms  1 . 

I 

IDG  j  Facility  radio  horizon  angle  in  degrees. 

IMN  |  Minutes 
1SEC  .  and  seconds 

DC I  j  Diameter  of  faclllcy  counterpoise  (ft). 

NCTE :  Zero  or  negsclve  values  will  cause  the  program 
i  to  assume  thac  no  counterpoise  la  present. 

HCI  Height  of  facility  counterpoise  above  facility  site  surft.ce. 

ICC  Code  for  counterpoise  reflection  material  type 
|  (Same  as  for  KSC  on  card  6.) 
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TABLE  B-2 


PARAMETER  FIELDS,  ATOA  (cont.) 


IIITRE  VERSIONS  1.243 
CARO  5 


Card 

Columns 

_  Symbol 

Description 

1-6 

HAI 

Height  of  aircraft  or  higher  antenna  above  mal 

7-9 

IAA 

Code  for  aircraft  antenna  pattern: 

(1)  Isotropic 

(2)  WF. 

(3)  TACAN  (RTA-2) 

(4)  4-loop  array  (cosine  vertical  pattern) 

(5)  8-loop  array  (cosine  vertical  pattern) 

(6)  1  or  11  (cosine  vertical  pattern) 

(7)  JTAC  with  tilt 

(8-21)  Special  antennas 

10-11 

JS 

Code  for  antenna:  (0)  directive,  (>0)  tracking 

12-13 

NPL 

Code  for  polarization  of  aircraft  antenna 

14-18 

T2T 

Tilt  of  the  aircraft  antenna  main  beam  In  deg. 

Not  used  for  patterns  1-6. 

19-23 

H2PBW 

Half  of  the  half-power-beamwldth  of  the  aircraft  antenna. 

Not  used  for  patterns  1-6. 

24-27 

ENO 

Surface  refractlvlty  referred  to  sea  level  (N-units)  from 
fit.  3.  NOTE:  301  N-unlts  will  be  used  If  value  Is 
specified  or  Is  <250  or  >400  N-unlts. 

28-33 

F 

Frequency  (MHz) 

34-39 

eirp 

Equivalent  Isotropically  radiated  power  (dBW)  for  power 
Density  and  power  available  output.  Sum  of  the  main 
beam  gains  In  dB  for  transmission  loss . 

40-45 

HPFI 

Elevation  of  effective  reflection  surface  above  mal. 
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TABLE  B-2.  PARAMETER  FIELDS,  ATOA  (cont.) 

MITRE  VERSIONS  1,  2  &  3 
CARD  6 


Card 

Columns 

Symbol 

1-2 

KSC 

- 

i 

j  TP 

8-12  I 

SCK 

13-14 

i 

ISS 

15-16  i 
| 

1  JM 

17-18 

1 

IOS 

1 

19-20  1 

I PK 

1 

1  1  _  *>0 

6  »  *4 

JO 

_ Description _ 

Code  for  earth  reflection  material  type  (table  2): 

(1)  sea  water 

(2)  good  ground 

(3)  average  ground 

(4)  poor  ground 

(5)  fresh  water 

(6)  concrete 

(7)  metallic 


Temperature  In  Celsius  (0°,  10° ,  20°)  of  the  water 

Sigma  In  ft  or  meters  If  you  do  not  wish  to  use  the  standard 
ones. 

Code  for  sea  state  per  table 

Code  for  sigma:  (0)  use  Standard,  (>1)  read  In  SCK 

Code  for  Ionospheric  Scintillation  Index  group*. 

j  (0-5)  See  Figure  5,  Coes  Report 

i  (-1)  Variable  group 

I 

|  Code  for  frequency  scaling  factor: 

I  (0)  not  use 

'  (1)  ( 136 / f ) 

! Code : 


(0)  no  scintillation 
(>0)  scintillation 


23-25 


KLM 


Code  for  climate  1  and  time  blocks:  (0)  Continental  all  years, 

(1)  Equatorial,  (2)  Continental  subtropical,  (3)  Maritime 
subtropical,  (4)  Desert,  (5)  Continental  temperate,  (6)  Mari¬ 
time  temperate  overland,  (7)  Maritime  tesiperete  overaeae, 

(8)  none,  (9)  Summer  time  block,  (10)  Winter  time  block, 
(11-18)  time  blocks  1-8,  (19)  All  year  time  block. 

(Code  (0)  is  recommended  for  best  results) 


26-28 


MX1  Code  for  mixing:  (0)  no  mixing.  Ary  other  number  is  weighting 
factor  for  climate  1. 


29-31 


KLM2  Code  for  climate  2  using  the  same  code*  as  those  in  KLM. 


32-14 


MX2  The  weighting  factor  for  climate  2  when  mixing. 


t 

i  i 


I 
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TABLE  B-2 .  PARAMETER  FIELDS,  ATOA  (cont.) 


.HllRE  VERSIONS  2  AND  3  ONLY 
CARD  7 


Card  I 
Co  lvijinsj _ I 


1-3  I 


4-6 


Qv  nb 
11 
12 


I 

H" - 


Description 


Code  for  percent  confidence 
Levels 


7-9  i 

13 

1 

;  (1)  0.0012 

(10) 

1.02 

(19) 

60.02 

(28) 

10-12  : 

14 

1  o  0.0022 

(11) 

2.02 

(20) 

70.02 

(29) 

i3-i5  ; 

15 

;  (3)  0.0052 

(12) 

5.02 

(21) 

R0.02 

(30) 

i 

16-18  , 

16 

.  (4)  0.0102 

(13) 

10.02 

(22) 

85.02 

(31) 

19-21 

’7 

!  (5)  0.0202 

(14) 

15.02 

(23) 

90.02 

(32) 

22-24  : 

18 

j  (6)  0.0502 

(15) 

20.02 

(24) 

95.02 

(33) 

25-27  , 

19 

!  (7)  0.10C2 

(16) 

30.02 

(25) 

98.02 

(34) 

28-30  1 

no 

i  (8)  0.2002 

(17) 

40.02 

(26) 

99.02 

(35) 

* 

n-33  | 

1 

in 

j  (9)  0.5002 

(18) 

50.02 

(27) 

99.52 

34-35 

ID 

i 

!  Code  for  heading 

99. 82 

99.92 

99.9S2 

99.987. 

99.992 

99.9952 

99.998% 

99.9992 


(0)  Will  produce  no  heading. 

(21  Heading  will  appear  at  the  beginning  of  every  page  of 
numerical  data.  The  heading  will  include,  title,  P«8' 
IZlll  date,  time,  and  partial  list  of  input  parameters. 
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TABLE  B-2.  PARAMETER  FIELDS,  ATOA  (cont.) 


ilITRE  VERSIONS  2  AND  3  ONLY 
CARD  S 


Card 

Columns 

Symbol 

Description 

1-4 

XCON  (1) 

First  line  of  sight  Increment 

5-8 

XCON  (2) 

Second  line  of  sight  increment 

9-12 

XCON  (3) 

Third  line  of  sight  increment 

13-16 

XCON  (4) 

Fourth  line  of  sight  increment 

17-19 

NTM  (1) 

First  line  of  sight  repetition 

20-22  ' 

! 

NTH  (2) 

Second  line  of  sight  repetition 

23-25  ' 

NTM  (3) 

Third  line  of  sight  repetition 

26-28  | 

NTH  (4) 

|  Fourth  line  of  sight  repetition 

29-32  1  YCON  (2)  | 

First  over  the  horizon  increment 

33-36  ; 

i 

YCON  (2)  | 

Second  over  the  horizon  lncremenet 

37-40 

| 

YCON  (3) 

Third  over  the  horizon  increment 

41-44 

YCON  (4)  | 

Fourth  over  the  horizon  increment 

45-47 

MTM  (l) 

First  over  the  horizon  repetition 

48-50 

MTM  (2) 

Second  over  the  horizon  repetition 

51-53 

MTM  ( 3)  1 

Third  over  the  horizon  repetition 

54-56 

MTM  (4)  j 

Fourth  over  the  horizon  repetition 

TABLE  B-3.  IMPORTANT  FILES 


ATOA  MITRE  VERSIONS 


ATOA.CLIST 
ATOAH3.CUST 
DEBE.CNTL 
FORT. CLIST 
LATOA. CLIST 
LATOAM3 . CLIST 
LISTOF? . CLIST 
MY AGAIN. FORT 
sample  3.  ern. 
SAMPLE. DATA 
SAMP  LEM.  DAT  A 
TAPE  .CNTL 


The  following  pages  contain  a  description  of  each  file. 
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TABLE  B-3.  IMPORTANT  FILES,  ATOA  MITRE  VERSIONS  (cont.) 


ATOA.CLIST 

ATOA.CLIST  Is  a  command  Use  that  executes  ATOA  MITRE  VERSION  1. 

To  use  this  cosmand  list  enter. 

EXEC  ATOA  'dnaae' 

dname  is  the  name  of  the  file  which  contains  the  lines  (cards) 
of  data  used  by  ATOA  MITRE  VERSION  L. 

The  oueput  la  stored  In  a  data  file.  The  name  of  the  output  file  will  be 
tne  name  of  the  input  file  with  a  "D"  suffix  attached. 

Example: 

SAMPLE. DATA  contains  card  set  up. 

You  enter:  EXEC  ATOA  'SAMPLE' 

The  output  Is  stored  In: 

SAMPLED. DATA 

This  method  of  executing  ATOA  Is  not  recomended  If  many  runs  are  to  be 
done.  The  preferred  method  is  submitting  the  runs  as  s  batch  job  (sea 
SAMPLE. CNTL). 

AT0AM3 . CL 1ST 

AT0AM3.CLIST  Is  identical  to  the  ATOA.CLIST  with  the  exception  that 
ATOA  MITRE  VERSION  3  la  executed  Instead  of  version  1. 

DEBE.CNTL 

This  file  la  a  tape  dump  control  file. 


TABLE  B-3.  IMPORTANT  FILES,  ATOA  MITRE  VERSIONS  (cont.) 

FORT.CLIST 

FORT.CLIST  execute*  the  Fortran  H  extended  compiler  by  utlllelng 
the  elgnif  leant  compiler  option*  AltTODBL(DBLFAD) ,  N0OPTIH1ZE,  and  ALC. 

ATOA  and  all  Its  subroutines  were  compiled  using  FORT.CLIST.  To  u»e 
thla  ccaund  list  enter: 

EXEC  FORT  'dname' 

dnaae  1*  the  nane  of  a  file  containing  Fortran  source  code. 

The  output  will  be  stored  In  an  object  file  under  the  same  name  as  the 
Fortran  source  file.  In  addition  a  taxt  fil*  will  b*  created  containing 
the  source  listing  with  ISN  numbers,  cross  reference,  and  diagnostic 
messages. 

Example : 

ATOA. FORT  contains  source  cod*. 

You  enter: 

EXEC  FORT  'ATOA' 

Theee  files  are  now  created: 

ATOA. OBJ  object  code 
ATOA. TEXT  diagnostics 

LATOA.CLlST 

LATOA.CH ST  Invokes  the  Linkage  Editor  and  turns  ATOA  MITRE  VERSION  1 
Into  an  executable  load  module. 

To  use  this  comund  list  enter: 

EXEC  LATOA 
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TABLE  B-3.  IMPORTANT  FILES,  ATOA  MITRE  VERSIONS  (cont.) 


the  load  nodule  la  stored  In  a  partitioned  data  aat  by  the  name  of 
ATOA.LOAD  (MODI) 

LAT0AM3.CLIST 

LATOAM3.CLIST  Invokes  the  Linkage  Editor  and  turna  ATOA  MITRE  VERSION  3 
Into  an  executable  load  nodule. 

To  use  thla  command  list  enter: 

EXEC  LATOAK3 

The  load  module  la  stored  In  a  partitioned  data  set  by  the  name  of 
ATOA.LOAD  (HOD3) 

L1ST0FF.CUST 

USTOFF.CLIST  will  create  a  complete  hard  copy  Hating  of  ATOA 
MITRE  VERSION  2. 

To  use  LISTOFF  enter: 

EXEC  LISTOFF 

MTACAIN. FORT 

MYAGAIN.FORT  is  a  complete  listing  of  the  original  ATOA  program  and 
all  Its  subroutines. 
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TABLE  B-3.  IMPORTANT  PILES,  ATOA  MITRE  VERSIONS  (cont.) 


SAMPLES. CNTl 

SAMPLES. CNTL  la  a  control  flit  chat  contains  all  the  necessary 
job  control  language  (JCL)  cosmndt  to  run  ATOA  KITRE  VERSION  3  plua 
ten  tinea  (cards)  of  data  uead  for  the  aaaple  teat  caae  of  ATOA. 

Thle  control  file  la  dealgned  to  aatva  aa  a  guide  to  prograonera  vho 
use  the  MITRE  version  of  ATOA. 

The  first  line  of  thla  file  la  the  JOB  statement  which  contains 
Information  that  la  needed  by  the  system.  A  complete  understanding  of 
the  JOB  statement  should  be  known  by  the  Prograasner  and  It  can  be  found 
in  the  Bedford  Computer  Center  Facility  Manual  on  page  2-5.  A  basic 
explanation  of  JCL  use  can  also  be  found  in  IBM's  FORTRAN  IV  Prograomer's 
Culde. 

The  second  line  la  a  continuation  of  the  first  line. 

The  third  line  of  the  control  file  contains  a  comment  which  la 
denoted  by  "//*"  in  the  first  three  columns.  If  one  aleeh  Is  removed 
(/ *  ROUTE  Pr.INT  HOLD)  the  Route  command  will  be  executed. 

The  route  conoand  allows  Job  output  to  be  viewed  at  the  terminal  before 
being  submitted  to  the  line  printer.  After  a  Job  Using  Route  is  completed 
the  three  most  used  time  sharing  option  (TSO)  commands  are: 

l)  L1STJES  Job-name 

Thla  command  displays  job  output  at  tha  tarminal.  Aftar  thle 
cosmuud  la  entered  many  subcommands  era  available  to  manipulate 
the  data.  These  comends  Include  PP,  PB,  FT,  fB,  END,  QUIT, 
DELETE.  A  complete  description  of  tha  commands  available  can 
be  obtained  by  entering  HELP  or  by  referring  to  tha  TSO/8UPERSET 
UTILITIES  user’s  guide. 


TABLE  B-3.  IMPORTANT  FILES,  ATOA  MITRE  VERSIONS  (cont.) 


2)  OUTPUT  Job-name  NOHOLO 

This  TSO  command  submits  a  Job  output  to  the  line  printer  and  Is 
generally  used  after  the  LISTJES  command. 

3)  CANCEL  Job-name  PURGE 

This  command  deletes  output  held  by  the  LISTJES  conoand  and  It 
Is  generally  used  If  the  OUTPUT  comma™!  ia  not  used  to  release 
the  Job  output. 

If  "//*"  Is  left  it.  the  first  three  columns  the  route  command  will 
not  be  executed  and  all  Job  output  will  go  directly  to  the  line  printer. 

The  fourth  line  of  the  control  file  gives  the  computer  the  name  of  the 
file  where  the  program  ATOA  can  be  found. 

The  fifth  line  Is  the  member  of  that  particular  file  which  is  to  be 
executed  by  the  computer.  H0D3  In  this  case  refers  to  ATOA  MITRE  VERSION  3. 

The  sixth  and  seventh  lines  allocate  input  and  output  for  the  Read 
and  Write  statements. 

The  eighth  line  sends  diagnostics  to  the  printer. 

The  ninth  line,  //3YSIN  DD*,  tells  the  computer  that  the  following 
lines  contain  the  Input  (card  set)  used  by  the  program.  In  the  particular 
case  of  ATOA  MITRE  VERSION  3  which  this  control  file  executes,  the  program 
expects  eight  lines  (card)  for  each  run.  Each  additional  run  of  the  program 
Is  done  by  adding  eight  more  lines  (cards)  directly  following  the  format  lines 
(cards)  used  for  the  previous  pass  of  the  program.  When  all  lines  (csrda) 
have  been  entered  in  multiples  of  eight,  two  additional  blank  lines  (cards) 
must  immediately  follow  for  proper  termination  of  the  computer  program  ATOA 
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TABLE  B-3.  IMPORTANT  FILES,  ATOA  MITRE  VERSIONS  (cont.) 


MITRE  VERSIONS  1,  2  &  3.  If  the  control  file  is  used  to  execute  version  1 
instead  of  version  2  or  3  the  only  difference  would  be  that  the  lines 
(cards)  would  be  entered  in  multiples  of  six  instead  of  eight  because 
of  the  fact  that  only  six  cards  are  required  for  each  pass  of  version  1. 

The  line  that  Immediately  follows  the  two  blank  lines,  a  slash  star 
(/*)  followed  by  blanks  is  standard  JCL  which  signifies  the  end  of  a 
data  block  (card  setup)  to  the  computer. 

The  very  last  line,  with  two  slashes  (//)  followed  by  blanks, 
signifies  the  end  of  a  Job. 

SAMPLE -DATA 

This  data  file  contains  th«  lines  of  data  (card  setup)  used  to 
run  the  ATOA  MITRE  VERSION  1  sample  case. 

SAMPLEM.DATA 

This  data  file  contains  the  lines  of  data  (card  setup)  used  to 
run  the  ATO«  MITRE  VERSIONS  2  and  3  sample  caeee. 

TAPE.CNTL 

This  control  file  le  used  to  enter  the  ATOA  program  from  tape 
Into  the  IBM  computer. 


